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We consider Laplacian in a straight planar strip with Dirichlet boundary 

' which has two Neumann "windows" of the same length the centers of which 

CN ■ 

are 21 apart, and study the asymptotic behaviour of the discrete spectrum as 
, Z — > oo. It is shown that there are pairs of eigenvalues around each isolated 

' eigenvalue of a single-window strip and their distances vanish exponentially 



in the limit / ^ oo. We derive an asymptotic expansion also in the case 
where a single window gives rise to a threshold resonance which the presence 
of the other window turns into a single isolated eigenvalue. 



> ■ 1 Introduction 

X 

\ Geometrically induced bound states in waveguide systems have attracted a lot of 

attention recently. The main reason is that they represent an interesting physical 
effect with important applications in nanophysical devices, but also in flat electro- 
magnetic waveguides - cf. |LCMj and references therein. At the same time, such 
a discrete spectrum poses many interesting mathematical questions. 

One of the simplest systems of this kind is a straight hard-wall strip in the plane 
with a "window" or several "windows" in its boundary modeled by switching the 
Dirichlet boundary condition to Neumann in the Laplace operator which will be 
the Hamiltonian of our system. By an easy symmetry argument it represents the 
nontrivial part of the problem for a pair of adjacent parall el waveg uides coupled 



by a window or several windows in the common boundary (ESTV ; this explains 
the name we use for the Neumann segments. 
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The discrete spectrum of such a system is nonempty once a Neumann window 
is present. Various properties of these bound states were analyzed including their 
number and behaviour with respect to parameters. Recently we discussed the 
way in which the eigenvalues emerge from the continuous spectrum as the window 
width is increasing .BEGl; refer to this paper for references to earlier work. 
Here we address a different question: we consider a strip with a pair of identical 
Neumann windows at the same side of the boundary and ask about the behaviour 
of the discrete spectrum as the distance between them grows. 

There is a natural analogy with the multiple-well problem in the usual Schrodin- 
ger operator theory - see [BCDj and references therein or |Da| Sec. 8.6] - even if the 
nature of the effect is different. Recall that in waveguides of the considered type 
there are no classically closed trajectories apart of the trivial set of measure zero, 
and likewise, there are no classically forbidden regions. Hence the semiclassical 
analysis does not apply here, in particular, there is no Agmon metric to gauge the 
distance of the windows which replace potential wells in our situation. 

Nevertheless, the picture we obtain is similar to double-well Schrodinger opera- 
tors. If the half-distance / between the windows is large, there is pair of eigenvalues, 
above and below each isolated eigenvalue of the corresponding single-window strip. 
We will derive an asymptotic expansion which shows that the pair splitting van- 
ishes exponentially as / — * oo together with the appropriate expansion for the 
eigenfunctions. On the other hand, the analogy a double-well Schrodinger opera- 
tor can be misleading. This is illustrated by the case when the single- window strip 
has a threshold resonance, which turns into a (single) isolated eigenvalue under 
influence of the other window. We derive the asymptotic expansion as / — >■ cx) for 
this case too; it appears that it is exponential again with the power determined by 
the term coming from the second transverse mode present in the expansion of the 
resonance wavefunction. 

Let us describe briefly the contents of the paper. In the next section we formu- 
late the problem precisely and state two theorems which express our main results. 
In Section 3 we collect general properties of the involved operator. Before coming 
to the proper proofs, we analyze in Sections 4 and 5 strip with a single window, 
in particular, we show how the original question stated in PDE terms can be re- 
formulated as a pair Fredholm problem, the second being obtained from the first 
one as a perturbation. Finally, in Sections 6 and 7 we prove Thms. 12.11 and 12.21 

2 Formulation of the problem and the 
main results 

Let X = {xi,X2) be Cartesian coordinates and suppose that H is a horizontal 
strip of a width d, i.e. H := {x : < X2 < d}. In the lower boundary of 
the strip we select two segments of the same length 2a. The distance between 
these segments, denoted as 21, will be large playing the role of parameter in our 
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Figure 1: Waveguide with two Neumann segments 



asymptotic expansions. We will employ the symbol 7/ (a) for the union of these 
segments, 7; (a) := 7/'(a) U 7f (a), where 7;"^ (a) = {x : |xi =F ^| < a, X2 = 0}. The 
remaining part of the boundary of 11 will be indicated by Ti[a) (cf. Figure P). 
The main object of our interest are discrete eigenvalues of the Laplacian in 11 with 
Dirichlet boundary condition on Ti[a) and Neumann one on 7^. We denote such 
an operator by Hi{a) and look what happens if / —>■ 00. 

In order to formulate main results of this paper we need some additional no- 
tations and preliminary results concerning a single-window strip. Denote 7(a) := 
{x : |xi| < a, X2 = 0}, where T{a) := dU \ 7(a). It was proven in (ESTVj that the 
Laplacian in 11 with Dirichlet condition on T{a) and Neumann one on 7(a) has 
(simple) eigenvalues below the threshold of the continuous spectrum for any a > 0; 
their number is finite and depends on a. We will indicate the operator in question 
and its eigenvalues by H{a) and A-,- (a), j = 1, . . . ,n, respectively, with the natural 
ordering, Ai(a) < A2(a) < ... < A„(a) < supposing that the corresponding 



eigenf unctions ipj are normalized in L^(n). Furthermore, it was shown in ESTV 
that there are critical values of size of Neumann segment, = oq < ai < 02 < • • • , 
for which the system has in addition a threshold resonance, i.e. the equation 
{H{an) + = has a nontrivial solution ip'^{x) unique up to a multiplicative 
constant. This solution and eigenfunction ipj mentioned above have a definite 
parity with respect to Xi and behave in the limit Xi +00 as 



^ljj{x) = a, e-/?-^"^ sin (^^) + O (^e 



(2.2) 



with some constants aj, it is clear that aj = aj{a). While normalization of 
ijjj is natural, the normalization of tp^ can be arbitrary, of course; we choose it in 
such a way that asymptotically the function coincides with the first normalized 
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transverse mode. Needless to say, when the window size is made larger than the 
critical value, the threshold resonance turns into a true eigenvalue. 
Now we are ready to formulate the main results. 

Theorem 2.1. Let the window length be non-critical, i.e. a G (a„_i,a„) for some 
n G N. Then the operator Hi{a) has for any I large enough exactly 2n eigenvalues 
\^(l,a), j = l,...,n, situated in the interval (^,^). Each of them is simple 
and has the asymptotic expansions 



(2.3) 



as I ^ oo for j = 1, . . . ,n, where a is an arbitrary fixed positive number. The 
coefficient fij is given by 



IJ,j{a) := aj{a) d\j — 



(2.4) 



or alternatively by 
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(2.5) 



The eigenfunctions ijj{x) associated with eigenvalues Xj{l,a), j = 1, . . . ,n, have 
a definite parity being even for A^(/,a) and odd for A~(/,a). Furthermore, in the 
halfstrips U"^ := {x : ±xi > 0, < X2 < d} they can be approximated by 



W^{U^) as I ^00. 



Theorem 2.2. Let the Neumann segment have a critical size, a = a^^. Then the 

2 2 

operator Hi{a) has 2n + l eigenvalues in (j^,^) for I large enough. The first 2n of 
them together with the associated eigenfunctions behave according to Theorem \2.1\. 
while the last one, X^,i{l,an), exhibits the asymptotics 



A^+i(/, a„) 
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where 



/i := 3/3^^2 , 



(2.6) 



(2.7) 
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or alternatively, 
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(2.8) 



r/ie associated eig enf unction ip^^-^ is even w.r.t. xi and for any R in the rectangles 
{x : \xi =F ^1 < R} n n can be approximated for large values of I as 



i:^^,{x) = r{x,Tl,x,) + 0{e-'^^) 



in Wo -norm. In addition it behaves in the limits xi ±oo as 



d \ 



X 
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Before proceeding further, let us recall what we have said in the introduction about 
the analogy with the multi-well problem for Schrodinger operators. As in that case 
a (simple) eigenvalue of the single window problem gives rise to a pair of eigenvalues 
(corresponding to eigenfunctions of different parities) which are exponentially close 
to each other with respect to the window distance, and moreover, in the generic 
case the splitting is determined by the eigenvalue distance from the threshold. 
At a glance the multiplicity is doubled by the perturbation, however, in reality 
the problem decomposes due to mirror symmetry into a pair of problems with 
definite parities whose eigenvalues tend to the same limit (see below and Sec. 5). 
On the other hand, the asymptotics ()2.fij) in the critical case differs from what 
the Schrodinger operator analogy would suggest being determined by the distance 
from the second transverse eigenvalue. 

Let us now describe our way to prove Theorems 12.11 and 12.21 The main idea is 
to reduce the eigenvalue problem at hand to a Fredholm operator equation of the 
second kind with a regular perturbation. Investigating this problem we will get the 
result both in the generic situation described in Theorem 12. II and for perturbation 
of a threshold resonance, just the analysis in the latter case is more subtle. 

Our task can be simplified by taking into account the symmetry of the problem 
with respect to reflections, xi —xi, which means that the operator decomposes 
into orthogonal sum of parts of a definite parity which can be considered separately. 
This allows us to cut the strip 11 into a pair of halfstrips 11^ and to consider the 
Laplacian in 11+ with Dirichlet condition everywhere at the horizontal boundaries 
of the halfstrip except for 7/^(0), where the boundary condition is Neumann. Ac- 
cording to the chosen parity of an eigenfunction we impose at that Dirichlet 
condition for odd eigenfunctions of the original problem at the vertical part of the 
boundary, xi = 0, or Neumann for the even ones. Moreover, it is convenient to 
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shift the halfstrip by Xi ^ Xi—l in order to fix position of the Neumann segment of 
the boundary. As a result, we arrive at the following pair of eigenvalues problems, 

-Aij = \ij, X G n' , 

dip (2.9) 
'0 = 0, X G r(a) , — — = , X G 7(a) , hu = , Xi = —I . 

0x2 

Here II' := {x G 11 : Xi > — /} is the shifted halfstrip and h is the boundary 
operator which acts as hu = u or hu = ^ in the odd and even case, respectively. 
Eigenvalues of ()2.9|1 obviously coincide with those of Hi{a) and by the even/odd 
extension one gets the eigenfunctions of the original problem. 

Finally, we remark that the problem has a simple behaviour with respect to 
scaling transformations which allows us to perform the proofs for d = it only. 

3 Preliminaries 

Let us collect first some general properties of the spectrum of our operators. 

Proposition 3.1. The discrete spectrum of the operator Hi{a) is non-empty for 
any I > a > 0. It consists of a finite number of simple eigenvalues contained in 
the interval (|, l) for d = ir which depend continuously on I and a; for a fixed a 
those corresponding to even and odd eigenfunctions are increasing and decreasing, 
respectively, as functions of the window separation parameter I. All eigenvalues of 
Hi{a) which remain separated from the continuum converge to those of H{a) as 
I +00, and to each eigenvalue of H{a) there exists a pair of eigenvalues of Hi[a) 
associated with eigenfunctions of opposite parities converging to that eigenvalue of 
Hi{a). If the Neumann segment has a critical width, a = an, then there is a 
unique eigenvalue (corresponding to an even eigenf unction) which tends to one as 
I — > +00. 

Proof. By the minimax principle and an elementary bracketing estimate the 
eigenvalues of Hi{a) can be squeezed between those of H{1 + a) and H{a). The 
essential spectrum of all the three operators is the s ame being equal to [1,cxd); 



this fact in combination with the results of |ESTVj shows that crdisciHi(a)) is 



non-empty, finite, and contained in l). A similar bracketing argument shows 
that the eigenvalues A^(/,a) of the problem ()2.9j) with Neumann and Dirichlet 
boundary condition at Xi = — /, respectively, satisfy A^(/, a) < Xj{a) < Aj(/, a) for 
j = 1, . . . ,n, where the upper bound is replaced by one if the Dirichlet problem 
has less than j eigenvalues. In fact, bracketing implies also the stated monotonous 
behaviour with respect to /, i.e. 

A+(r,a) < A+(/,a) < Xj{a) < Xjil,a) < \]{l',a) (3.1) 

for V > I with the same convention as above; it is sufficient to write Uu as a union 
of and a rectangle separated by an additional Neumann or Dirichlet boundary 
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condition and to realize that in neither of these cases the rectangle can contribute 
to the spectrum below the continuum threshold, because it has Dirichlet condition 
at the horizontal part of the boundary. In addition, the standard domain-changing 
argument |0 Sec. VII. 6. 5] shows that the functions A^(-, a) are continuous. In view 
of the monotonicity mentioned above their limits as I —>■ oo exist; it remains to 
check that Xf{l,a) Aj(a). 

Take ip G D{H{a)) and a function g G such that g{x) = for x < and 
g{x) = 1 for a; > 1. Denoting hi{x,y) := g{2{x + /)//) we can construct a family 
{ipi} by ipi{x,y) '■= 4'{x,y)hi{x,y); by construction the function ipi belongs to the 
domain of H^{a) which is the Laplacian with the boundary condition as in ()2.9j] 
for hu = u. Using the fact that ||V/i;|p = 2||(y('|p/~^ and ||A/i;|p = 8||(7"|p/"^ one 
can check easily that ipi ^ ip and Hf{a)ipi H{a)ip as / ^ oo, so Hj^ {a) H{a) 
in the strong-graph sense. By |RSt Thm VIII. 26] this is equivalent to the strong 
resolvent convergence, hence to each Xj{a) there is a family of A~(a) converging 
to that value. Since the spectrum of Hi{a) in (i, 1) is discrete, simple, finite, and 
depends monotonously on /, we get the desired result. In a similar way one can 
check that A^(Z, a) Xj{a) as / ^ oo. 

The continuity w.r.t. a is proved as the a-continuity in case of a single Neumann 
window. We expand the solution inside and outside the window regions with 
respect to the appropriate transverse bases and match the Ansatze smoothly at the 
window edges. This yields an infinite family of linear equations for the coefficients 
of the expansions, which can be regarded as a search for the kernel of a certain 
operator in the i"^ space of the coefficients with a properly chosen weight. One has 
to check that this operator is Hilbert-Schmidt and continuous with respect to the 
parameters in the Hilbert-Schmidt norm. The argument is analogous to that from 
the proof of Proposition 2.1 of ^BEG . so we skip the details; the only difference 
is that due to the lack of symmetry the matching has to be performed at each 
window separately and the coefficient space is "twice as large" . 

It remains to check the last claim. Using bracketing once more we see that if 
the presence of the other window turns a threshold resonance into an eigenvalue, 
the corresponding eigenfunction must be symmetric; in view of the proved mono- 
tonicity it is sufficient to show that this happens for / large enough. Since this part 
of the proposition is not used in the proof of the claim of Theorem 12.21 concerning 
the first 2n eigenvalues A^, we may assume that it is already proven and that we 
thus know that for large / the operator Hi(an) posesses 2n eigenvalues A^ G (|, l) 
corresponding to eigenfunctions ip^. We seek a (2n-|- l)-dimensional subspace such 
that for any from it we have {ip, Hi{an)ip) — ||'0||^ < 0. To this aim we employ a 



Goldstone-Jaffe-type argument inspired by |ESTVj and choose 



where ip'^ is the resonance function 1)2.11) . xl,<; : M — > (0, 1] equals one in {—L,L) 
for some L > I + a and xl,? = exp(— <^(|x| — L)) otherwise, and p is a function 
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supported in the other window region. In view of the asymptotic behaviour ()2.ip 
such functions span a subspace of the needed dimension. Evaluating the energy 
form {ip, Hi{an)ip) — WipW^ we see that if some of the coefficients cf is nonzero, it is 
negative even with e = 0. In the opposite case we use the fact that in the leading 



term we have, as in [ESTVj . two competing terms, one linear in e and the other 
positive coming from the tails of controlled by the parameters L and we can 
choose them in such a way that the form is negative again. 

Hence Hi{an) has for / sufficiently large at least 2n + 1 eigenvalues. In fact, 
it has exactly this number, because its symmetric and antisymmetric parts have 
for large / enough / + 1 and / eigenvalues, respectively, otherwise we would have 
an contradiction with the monotonicity and continuity properties stated above. In 
particular, the largest eigenvalue is increasing w.r.t. I since it corresponds to an 
even eigenfunction. In view of ()3.1|) and the fact that a„ is the critical width, we 
conclude that A^_,_^ — 1— as / — oo. | 



4 Analysis of the limiting operator 

After these preliminaries let us pass to the proper subject of the paper. First we 
are going to discuss the limiting, i.e. one-window operator which means to analyze 
the following boundary value problem, 

Ou 

-(A + A)m = /, xen, M = 0, x G r(a) , — = 0,xe7(a). (4.1) 

0X2 

The right hand side / is here assumed to be finite and to belong to L^(n); our 
aim is to discuss the existence and uniqueness of the solution to ()4.1|) as well 
as its dependence on A. The method we use is to reduce M.lj) to a Fredholm 
operator equation. Then our task will be reduced to analysis of operator families, 
in particular their holomorphic dependence on the spectral parameter A (one need 
not specify at that the topology - cf. |R,S[ Sec. VI.3]). The reduction will follow a 
general scheme proposed by Sanchez-Palencia |SPj and it will be analogous to the 
treatment of a similar problem in |BEG| Sec. 3.1]. 



We will use the symbol Vs to indicate the open subset {A : ReA < 5} of the 
complex plane. The structure of the solution to the problem ()4.H1 for A close to 
1 and for A separated from 1 is different. This is the reason why we will consider 
these two cases separately. We suppose first that A e P^, where A„(a) < 5 < 1. In 
this situation it is sufficient to consider solutions of the problem ()4.H) in the class 
of functions which behave as 0{e~^^~''^^^^^) in the limit — 00. 

Since the function / is finite by assumption, its support lies inside the rectangle 
lib := n n {x : la^il < b} for some 6 > 0. Consider two boundary value problems, 

-{A + X)v^=g, xeU^, v^=0, xedll^, (4.2) 

where g is an arbitrary function from i^^(n) with the support contained in Ua 
for some A > max{a, 6 — 1}. This choice is given by the requirement that 11^ 
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contains both the window and the support of / in such a way which will make the 
smooth interpolation ()4.6|) used below possible. The problems ()4.2|) can be easily 
solved by separation of the variables; using the explicit form of Green's function 
of Laplace-Dirichlet problem on a halfline we get 



w± (x) = j G (x, t, A) g{t) dH , (4.3) 
n± 

oo ^ 

G^{x,t, A) = y -— _ gT'^,{A)(xi+ti)\ ^ ^4 



where Kj(A) := a/j^ — A . In the following we will also employ the "glued" function 

V equal to if a;i > and to v~ if xi < 0. The functions can be naturally 
regarded as results of action of the bounded linear operators T^{\), i.e. we have 

= T^{X)g, where : ^^(n^) iy|(n±) with the "halved" rectangles = 
n n {a; : < ±Xi < A}. It is easy to check that the operator families are 
holomorphic in A G P^. In the next step we consider the problem 

dw 

Aw = Av , xEUa, t: — = 0, a; G 7(a) , w = v, x G 911^ \ 7(a) . (4.5) 

0x2 

The function v may have according to its definition given above a weak discontinu- 
ity, i.e. a jump of the first derivatives. Thus we have to say what we mean by Av 
in ()4.5j) : it is the function from L'^(n.) which coincides with Af+ if Xi > and with 
Av~ if Xi < 0. With the problem ()4.2|) in mind we can also write Av = ~{Xv + g). 
The problem ()4.5|) is posed in a bounded domain, hence the standard theory of 
elliptic boundary value problems is applicable. In particular, we can infer us- 
ing |Laj that the function w exists, it is unique and belongs to iy2HnA)- We 
will also consider its restriction avoiding the points where the boundary condition 
changes, regarded as an element of VF2HnA) H 11^2^(11^ \ 5*^) for any r > 0, where 
Sr = {x : (xi ±a)^ + X2 < r^}. In this way we introduce a linear bounded operator 
T2 : L^illA) WliliA) n W^illA \ Sr) (for any r) such that w = T2C/. 

Next we employ a smooth interpolation. Let x be an infinitely different iable 
moUifier function such that %(''') = 1 if |t| < A — 1 while for |r| > Ait vanishes. We 
will construct a solution to the problem ()4.1|1 interpolating between the functions 

V and w, specifically 

m(x) = x{^i)'^{^) + (1 ~ x{^i))v{x). (4.6) 

Since w = T2g and = T^{X)g, we can also regard u as the result of an action 
of some linear operator T^i^X) which maps L'^(IIa) into 1^2^^^) H 14^2^(11 \ 5"^) for a 
fixed r > 0. Such an operator is linear and bounded, and as an operator family 
with respect to A it is again holomorphic. 

Owing to the definition of w and v the function u satisfies all the boundary 
conditions involved in 1)4.11) . and consequently, it represents a solution to 1)4.1)) if 
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and only if it satisfies tlie differential equation in question. Substituting ()4.6p into 
the latter and taking into account ()4.2|) . ()4.5|) . we arrive at the equation 



g + T,{X)g = f, 



(4.7) 



where T4 : L'^{JIa) L"^ 



(n^) is a linear bounded operator defined by 



T^{\)g := -2V,x ■ V,(^i; -v)-{w-v){A + \)x , 



(4.8) 



where the dot in the first term denotes the inner product in M^. The relation (j4.7|) 
is the sought Fredholm equation, considered in the space L'^{IIa)- Naturally the 
first thing to do here is to check the compactness of the operator T4. It can be done 
as follows. The function w — v belongs to 1^2^(11^), thus the operator mapping g 
into w — v is bounded as an operator from L'^(JIa) into W^2^(n^), and consequently, 
it is compact as an operator in the space L'^{IIa)] this solves the question for the 
second term at the right-hand side of ()4.8|) . Furthermore, due to the definition 
of the moUifier x the support of V^X lies within 11^ \ Ha-i- This domain does 
not contain the endpoints of the segment 7. Hence w — v & iy2^(supp V^^x); 
therefore Vx{w — v) considered as an element of L^(supp V^x) results from action 
of a compact operator mapping L^(n^) into L^(suppVxX)j this concludes the 
proof of compactness T^IX) considered as an operator in the space L^(n^). In a 
similar way one can check that T4(A) is a holomorphic operator family w.r.t. A. 

This conclusion allows us to apply to ()4.7|1 the standard Fredholm technique; 
we will see that solution to ()4.7|1 exists and is unique for almost all A except for 
points where a nontrivial solution for ()4.7p with zero right-hand side exists. This 
will yield a solution to our original problem because the two are equivalent; this 
is the contents of the following lemma the proof of which we skip because it is 
completely analogous to that of Proposition 3.2 in [BECilj . 

Lemma 4.1. To any solution g of there is a unique solution u = T3{X)g 

of \4.1\ ), and vice versa, for each solution of there exists a unique g solving 
\4- 7| ) such that u = Ts{X)g. The equivalence holds for any X ^Vs- 

Thus the equation ()4.7|) says how to find a bounded solution to ()4.ip : one should 
solve the equation ()4.8j) and then to construct the solution of ()4.1|1 by the procedure 
described above, i.e. by putting u = T3{X)g. 

Since the operator family T4(A) is holomorphic, the corresponding resolvent 
family (J + T4(A))~^ is meromorphic and its only poles are exactly the eigenvalues 
of H{a) - cf. |iSP^ Chap. 16, Th. 7.1]. In order to prove Theorem 12. H we need to 
know more about the behavior of (/ + T4(A))~^ in the vicinity of these poles. 

Lemma 4.2. Let Xq < 1 be an eigenvalue of H{a). Then for any X close enough 
to Ao the following representation is valid, 



(i + UX)) 



A- A, 



n + TeiX), 



(4.9) 



-0 
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where T^f := —{f,'ip)L^(u) andTQ : L^(n^) L^(n^) is a hounded linear operator 
which is holomorphic in A. Furthermore, is such that ip = T3(Ao)0, where ip is 
an eigenfunction of H{a) associated with Aq and normalized in L^(n). 

Proof. We assume throughout that A G 1^5 hes in a small neighborhood of Aq 
containing no other eigenvalues of H{a). As we have already mentioned, the 
operator family (/ + T4(A))~^ has a pole at Aq. It means that the vector-valued 
function g : X {I + T4(A))~^/ satisfies 



9~q 

(A - Ao)^ (A - Ao)-?- 



^(A) = 7V^ + 7V^^, (4.10) 



where g is a positive integer and g is holomorphic in A. Substituting this represen- 
tation into ()4.7|1 and calculating the coefficients of (A — Aq)"'' we see that g-g must 
satisfy the equation g_g + T{XQ)g_g = 0, in other words g^g = (pT^f, where T^f is 
a number depending on /. Together with ()4.10j) this means that the solution to 
()4.1|) associated with g, i.e. u = T2,{\)g, can be written as 

where is u is holomorphic in A. Due to the definition of T^IX) this formula is valid 
in the sense of iy2^(n)-norm as well as in W^2^(n \ Ha)- Taking the inner product 
of ()4.ip with ip, using the fact that the latter is an eigenfunction of H{a), and 
performing an integration by parts in 11^ with R large enough we find 



The functions u and ip behave at infinity as (9(e~'^^'^'^^^) and (9(e^'^'^'^'^^^^), 
respectively. With this fact in mind we can pass to the limit i? — ^ oo in ()4.12|1 for 
each fixed value of A; this implies the identity 

Ao(M,V^)L2(n) = (/,V^)L2(n) + A(M,V^)L2(n) . 

Substituting to it from ()4.11|1 and computing the coefficients at the same powers 
of A — Aq, we see first that g = 1, and furthermore, that T5/ = —{f^ip)i,2(jx)- This 
completes the proof | 

We will also need to know the behavior of the inverse (/ + T/^{X))^^ as A ^ 1. 
For the right hand side / in 1)4.11) with a definite parity w.r.t. xi = it was done 
BEGj . here we have just to show how to extend this result to our case. We 



m 



will assume that A lies in a small neighborhood of one and that this neighborhood 
contains no eigenvalues of H{a). First of all, however, we should characterize the 
class of functions in which we will seek the solution of 1)4.1)) in this case. Instead 
of A we introduce another parameter by setting A = 1 — k^, where k lies in a 
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small neighborhood of zero. The only restriction to the size of this neighborhood 
correction is that the associated values of A should not coincide with eigenvalues 
of the operator H{a). If k is real and a solution to the problem ()4.H) exists, it is 
unique and holomorphic in k. This fact follows from the arguments given above, 
because for such a A = 1 — the equation ()4.7|1 is uniquely solvable. The said 
solution can be extended to all values of k in the vicinity of zero so that this 
extension will be an analytic function of k. The existence of such an extension is 
guaranteed by the definition of the functions in ()4.3|) where Ki{X) is nothing 
else than k introduced above. We see that the formulae ()4.3|1 are valid not only 
for real k but also in a complex neighborhood including k = 0, because the kernels 
()4.4j) have finite limits = 0, namely 

G^{x, t, 1) = (|xi — til =F (a^i + ti)) smx2 sint2 

1 (4.13) 

+ V — (e-'^^(m^i-*il - eT-.(i){^i+ti)) sin jX2 sin jt2 . 

This is why we are able to extend the solution of the problem ()4.H) analytically 
to all values of k in the vicinity of zero. We should also stress that the function 
u given by ()4.6|) decays exponentially at infinity if Re k > 0, it is bounded for 
Re/« = and increases exponentially provided Re/t < 0. 

In this approach all the operators introduced above preserve their properties 
when we vary the range of the variables passing from unbounded domains the cut- 
off ones treating, for instance, as operators mapping L^(n^) into 1^2^(11^) for 
any R. As another notational simplification we will not introduce an extra symbol 
for the composed mapping k, ^-^ Tj(l — k^) and write instead just Ti{K). 

Mimicking the argument used in the proof of |BE(t| Thm 3.4], one can check 
the following claim: 

Lemma 4.3. // the Neumann segment of the boundary does not have a critical 
size, the operator {I + Ti{K))'^ exists and is uniformly hounded in n in the vicinity 
of zero. In the opposite case, i.e. a = an, we have in a punctured neighborhood of 
zero the following representation, 

{I + T,{K)r' = ^T, + Ts{K), (4.14) 

K 

where Tjf := | (/,'?/'"' )L2(n) andT^ : L^(n^) — > L^(n^) is a bounded linear operator 
which is holomorphic in n. Furthermore, 0" is such that ip"' = T^{k = 0)0", where 
tp^ solves the equation {H{an) + = and behaves at infinity in accordance 
with 
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5 Analysis of perturbed operator 



The main purpose of this section is to reduce the problem 

-(A + A)m = /, xeU\ 

dip (5.1) 
ip = , xET(a), — — = 0, X G 7(a) , hu = , Xi = , 

ox I 

to an operator equation similar to ()4.7|1 . We will show that the problem 1)5.11) can 
be reduced to solution of a Fredholm equation which is a regular perturbation of 
the equation ()4.7|) . We will start from the case of Dirichlet condition at the cut 
Xi = 0, i.e., hu = u. We are going to employ the same scheme as in previous 
section and use the same notations unless stated otherwise. 

First we will treat the case X E Vs. In analogy with ()4.2j) we consider two 
problems, 

- (A + A)t;+ = ^ , X G n+ , = 0, X e (9n+ , (5.2) 

- (A + X)vi' = g, X eUY , vf = 0, X e dUY . (5.3) 

The first one coincides with ()4.2j) for v'^ , while in ()5.3j) we take into account the 
perturbation. Consequently, we have v^' := v~^, where is the function from ()4.3j) . 
The problem ()5.3)1 differs from ()4.2j) but it can be solved again by separation of 
variables. It is convenient to write its solution v'[' in the following form, 

v-{x)=v-{x) + j G^{x,t,\)g{t)dh, (5.4) 
n- 

G, (x, t. A) = — > -— -— — sinhK,(A)xi sinh K, (A)ti sin jx2 smjt2, 

(5.5) 

where v~ is given by ()4.3)1 : we keep in mind here that g is finite, and therefore its 
support lies inside 11; for all / large enough. As in previous section we can introduce 
a linear bounded operator Tq{X) : L^(n^) 1^2^(11;") such that = TQ{X,l)g. 
This operator can be represented as the sum Tg{X,l) = T^{X) + Tio(A, /), where 
Tio(A,Z) : L^(n^) — > 1^1(11;") is holomorphic in A, jointly continuous with respect 
to (A, /) provided A G Vs, I G [lo, +00], and Iq is a fixed number large enough. The 
norm of the operator Tiq is of order 0(e~'^^~'*') as / — > +00, hence for X E Vs we 
may consider this operator as an exponentially small perturbation. 

The analogue of the function w (denoted here by wi) is defined as above without 
any changes, i.e. as a solution of the problem ()4.5)1 with v replaced by 



vi :-- 



Vi , xi > , 
Vj^ , Xi < . 
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The solution of ()5.1|) is then constructed as an interpolation ()4.6|) with v and w 
replaced by vi and wi, this leads us to the desired operator equation, 



g + T^{\)g + T^^{\l) = f . (5.6) 

Here T4 is the operator appearing in ()4.7|) and Tii(A,/) : L^(n^) — > L^(n^) is a 
compact linear operator which is holomorphic in A and jointly continuous w.r.t. 
(A, /) provided A G and / G [Zq, +00]. The norm of the last named operator is 
exponentially small as I —>■ +00 uniformly in A G T>s: 

llTnll = 0(e-2'v^). (5.7) 

The solution to the problem 1)5.11) can be reconstructed from the function g by 
u = T3(A)(7+Ti2(A, l)g, where T12 : L^(n^) — > 1^2^(11') is a linear bounded operator 
the norm of which satisfies 

IIT12II = C»(e-'^^). (5.8) 

This operator is also holomorphic in A and jointly continuous with respect to 
(A,/) E Vs X [\q,+oc]. The equation ()5.(i|l is a second-kind Fredholm operator 
equation and it is equivalent to the problem (jOJ); this claim can be checked in 
the same way as we did it for ()4.7p in the previous section. 

The case of hu = ^ is treated in full analogy. The only difference due to 
another boundary condition at Xi = is the definition of the operator Tio which 
is now described by the kernel 

G, (x,t,X) = > --— ; — — - sinh /t,(A)a;i sinh/€,(A)ti sin 7x2 sin?'t2- 

^ nKj{X) cosh Kj{X)l 

(5.9) 

All the arguments used above remain valid. 

On the other hand, for A in the vicinity of one almost all the above arguments 
remain valid provided we replace A by (1 — /t^). In analogy with the previous 
section the operators introduced here may be considered on cut-off strips, i.e. as 
: L\U-^) ^ WiiU],), Tio(«:,0 : ^^(n^) ^ WHU],), TuinJ) : L2(n^) ^ 
L^(n^), Ti2(k, I) : L^(JIa) — > W2(JIr) for any fixed R. However, we are not longer 
allowed to say that these operators are holomorphic in k because of the terms 



sinh nl cosh kI 

in ()5.5j) . ()5.9|1 . since these terms have poles at k = yj and n = y(j + |)- Moreover, 
these terms are also responsible for the fact that the operators have no proper limit 
as K ^ and / +00. At the same time, restricting the range of k we will be 
able to show that the operators Tio, Tn, and T12 are small for small k and large /, 
thus we will be allowed to consider them as small perturbations again. This claim 
leans on the following lemma. 
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Lemma 5.1. Let k G (0, ^) he fixed and Q,^ := {k : | argfi;± || > x}. Then there 
is C > such that for small k G and large I the following estimate is valid, 



max 



K,e 



sinh kI 



K,e 



cosh kI 



<C{\k\+1-^) . 



Proof. We will show how to derive the first estimate, the proof of the second one 
is similar. We start by introducing the function 

P(z):=- ^ 



e^ - 1 

Suppose that z G Q^. If we have in addition \z\ < 1, one can check that 

\P(z)\<C 



(5.10) 



with some C independent on z. On the other hand, if > 1, z G Q^, and 
Rez > 0, then the exponent in the function P increases as \z\ oo and we arrive 
at ()5.1Up again (in general with another C). Finally, if |z| > 1, 2; G Q^, and 
Rez < then the exponent in the function P decreases and we have a uniform 
estimate, 

\P{z)\ < C\z\. 

Combining it with ()5.10p we get the inequality 

\P{z)\<Ci\z\+C2 
valid for 2 G Q and suitable Ci, C2. The obvious identity 

-kI 



P{2kI) 



sinh kI I 

then completes the proof of the lemma. | 

Using this result one can check that the operators Tio, Tn, and T12 are small 
for small k G Q>c and large /, holomorphic in k, and jointly continuous in (k, I). 



6 Proof of Theorem 12.11 

In this section we are going to derive the asymptotic expansions for the eigenvalues 
of Hi{a) separated from the continuum. We will also find the asymptotic behavior 
of the associated eigenfunctions. 

The main idea behind the calculation of the asymptotics is borrowed from 
|GaHIGa2llBEGj . Instead of dealing with eigenvectors of Hi{a) directly we consider 
here those of the problems ()2.9|) . In order to find eigenvalues of the latter we should 
look in accordance with the results of the previous sections for A such that the 
operator equation 

<l> + T4(A)<l> + Tn(A,/)$ = (6.1) 
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has a nontrivial solution. We will deal with eigenvalues which are close to a fixed 
eigenvalue Xj{a) of the limiting operator H{a); for simplicity we will denote the 
latter as Aq in the following. Also the parameter A will be assumed to be close to 
Ao, more specifically, it will be supposed to lie in a neighborhood of Aq containing 
neither any other limiting eigenvalue nor the point A = 1. 

By the definition of Tn the term Tii(A,/)$ in (jfi.lj) is supported inside 11^. 
Hence considering it as the right hand side, we arrive at the equation ()4.7|) with 
/ = — Tii(A, Choosing A 7^ Aq, we can invert the operator / + T4(A) obtaining 

$ + (/ + T4(A))-iTn(A,/)<l> = 0. 

Using the Lemma f4. 21 we can rewrite the last equation in the form 

$ - T^(^' Wlhu) + n{X)Tn{X, /)$ = ; (6.2) 

recall that G -^^^(11) here is the normalized eigenfunction associated with Aq and 
(f) e L'^(JIa) is a function such that ip = T3(Ao)0. 

The operator Tii(A, /) is small in the asymptotic region, / +00, while Tq{X) 
is holomorphic in A. Thus we may invert the operator / + Te(A)Tii(A, I) and apply 
the result to the equation ()6.2|) . which then acquires the form 

$ - Tn(A, mLHu){I + Te{X)Tn{X, /))" V = . (6.3) 

A — Ao 

The inner product (■?/', Tii(A, /)$)L2(n) does not vanish. Indeed, otherwise the func- 
tion $ would be zero too, however, we seek a nontrivial solution of the equation 
()6.H) . With this fact in mind, we express the function $ from the equation ()6.3p 
and then calculate the inner product (■?/', Tii(A, /)$)2,2(n). This procedure leads us 
to the equation 

^ - J^To ^^^^ + ^6(A)Tn(A, /))-V)^2(n) = , 

or in a more convenient form 

A - Ao - {^P, Tn(A, /)(/ + n{X)T^,{X, /))" V)i2(n) = • (6.4) 

This is the sought equation determining the perturbed eigenvalues of the problem 
()2.9j) . and, thus, of the operator Hi{a). The associated solution of the equation 
()6.H) . as it follows from ()6.3p . can be written as 

$ = (/ + T6(A)Tn(A,/))-V; (6.5) 

we naturally keep in mind the fact that the eigenf unctions are defined up to a 
multiplicative constant. 

The equation ()6.4j) determine all eigenvalues of Hi{a)] due to the equivalence 
between fl6.1|l and ()2.9j) only the eigenvalues of Hi{a) satisfy this equation. Thus, 
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by Proposition I3.H for every Tn there exists an unique solution of the equation 
()6.4p converging to Aq as / — +00. 

The desired asymptotic expansions for the perturbed eigenvalues can be calcu- 
lated directly from the equation ()6.4|1 . First of all we recall the assertion ()5.7|1 which 
implies that for A close to Aq the norm Tn can be estimated by (!?(e~*^^^^~'^°~'^^'). 
It allows us first to establish the estimate 

A - Ao = C (^e'i^VT^-^)i^ ^ 

and secondly to expand the second term in the equation ()(i.4|l obtaining 

A - Ao - iiJ, Tn(A, /)0)^.(n) + ^(e-^^^^--)') = . (6.7) 

We can also extract the leading term from the operator Tii(A, /), which obviously 
comes from the lowest-mode contribution to the sum at the right hand side of 
()5.5p . We will do that for hu = u, in the other case one proceeds analogously. 
First we introduce additional notations setting 

^g-2Ki(Ao)i I" 

I — -— sinh fti(Ao)a;i sin X2 / sinhfi;i(Ao)tisint20c^^^ , Xi < , 

V{x) := { ^f^ii^o) 



n- 

0, a;i>0. 

(6.8) 

Suppose that a function W solves the problem ()4.5j) with v = V, then 



Tn(A, /)0 = -(A + Ao) (V + xiW-V)) + 0(/e-2(2v^-)') in ^^(n 



A 



Using this identity together with the fact that the function Tii(A, /)0 is finite, we 
can calculate the leading term of the second summand in ()6.7|) . 

(^, Tn(A, /)0)i2(n) = - j {A + Xo) {V + x{W-V)) (fx + 0{le 



■2{2y/T^~a)l\ 



n 



{A + \^){y + x{w-v))(fx= liui [ (^p^-v^]ds 

n {x: |a:i|=_R, 0<a:2<7r} 

ds . 



-R^+^oo j dxi ^ dxi) 



{x: xi=—R, 0<X2<7r} 

(6.9) 

In order to calculate the last integral we use the fact that in view of the relation 
ip = T3(Ao)0 and the definition of T3 the constant a = aj in ()2.2|1 is given by 

2 f 

a = — -— / sinh(Ki(Ao)ti) sint2 0(i) c^^i , 

7i"fi:i(Ao) J 
n- 
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where p is 1 if even and -1 if it is odd. Using this relation together with 
and ()6.8|) . we can finish our calculations in ()6.9p arriving at 

lim / ( - V^] ds = na'K^iXo) e-^-iC^o)/ . 

R^+oo J \ dXi OXi J 

{x: xi=—R, 0<a;2<vr} 

Combining this with ()6.9|) and ()6.7|) we get the asymptotic ()2.3|) . ()2.4|) for Xj{a). In 
the case hu = a similar reasoning leads to asymptotics ()2.3|) . ()2.4|) for A^(a). 
In order to prove relation ()2.5j) it is sufficient to express a in terms of suitable 
integrals. Keeping the parity of ip in mind we compute 

0= lim e^'^^'''>'''smx2{A + Xo)ilj{x)d'^x 

R-^+oo J 

= lim f ( e'^i(^«)^isina;2^V(a;) - V'(a;)^e''i(^o)^isinx2 I cis 

R^+oo J \ OU OU J 

ij{x) e^^^^i dxi - an^l-Xo- 



7(a) 

This result leads us to formulae ()2.5p . 

The asymptotics of the eigenfunctions can be derived easily. The definite parity 
of those associated with X^{l,a) is obvious. The relation ()6.5|) tells us that 



$± = + C(e-('^^^-")). (6.10) 

The symbol "±" indicate here two variants of definition of the operator Tn. Now 
in order to prove the expansions for the eigenfunctions one has just to use this 
expression and to employ the arguments of the previous two section. More pre- 
cisely, we have ijj = T3(Ao)0 and = (T3(A^) + Ti2(A^,0)*^, where is the 
eigenfunction of the problem ()2.9p associated with the chosen eigenvalue and cho- 
sen variant of boundary operator h. Using ()fi.l()|l and the holomorphy of T3, T12, 
the estimates and ()5.8j) . we arrive at the asymptotical formula 

^± = + C fe-'(^^"")') (6.11) 



in 1^2^ (n'). Recovering now the eigenfunctions of Hi{a) we obtain all their prop- 
erties stated in Theorem 12.11 

Let us finally prove that there are no other eigenvalues of Hi{a) in Vi. Consider 
the equation ()6.1|) where A is close to one and does not lie in real semi-axis [1, +00), 
more specifically, suppose that k G Qj^. Then we can invert the operator (/ + 
T4(k))^^, and arrive at the equation 

$ + (/ + T4(/s:))-iTn(K,/)$ = 0, 
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where the operator (J + T^^k))'^ is uniformly bounded in k, because the Neu- 
mann segment does not have by assumption a critical size - see Lemma 14.31 
- while Tii(k, /) is small for all possible values k and /. Hence the operator 
(/ + T4(k))~^Tii(k, /) is also small, and therefore we can invert in turn the opera- 
tor (/+(/ + T4{k))~^Tii{k,, /)) which immediately leads us to the unique solution 
$ = 0. Moreover, the operator Hi{a) cannot have eigenvalues corresponding to k 
satisfying | arg/t ± ^| < x. Re k 7^ 0, simply because it is self-adjoint and all its 
eigenvalues are real, thus there is no other eigenvalues to Hi{a) in Vi. By this the 
proof of Theorem 12.11 is complete. | 

7 Proof of Theorem [212 

It is sufficient to consider in detail only the eigenvalue A„+i emerging from the 
continuum because all the statements related to the other eigenvalues verify in a 
way completely analogous to the previous section. 

We know from Proposition 13.11 that the eigenfunction associated with the in- 
dicated eigenvalue is even with respect to Xi, thus we have to consider here only 
the case hu = Assuming k E Q^^, we start with the equation 

$ + T4(«:)$ + Ti3(/€,/)$ = 0, (7.1) 

which is how ()6.1|) looks like in the present case, with Ti3(k, /) being the perturba- 
tion operator associated with ()5.9|1 . The operator Ti3(k, I) is small by Lemma EHl 
and an argument analogous to that which lead us to ()(j.3j) yields the equation 

$ + ^(^, Ti3(/€, mLHu){I + T8(A)Ti3(A, /))- V = . (7.2) 

Recall that = ip"" is a. solution to the equation (if(a„) + l)'?/'" = which behaves 
at infinity in accordance with 1)2.11) and (p ^ L'^O^a) such that = T^Ik = 0)0. 
From this equation one can deduce an analogue of the equation ()(j.4)l . namely 

2/€+ (^,Ti3(/s:,/)(/ + T8(K)Ti3(K,/))-V)i2(n) =0. (7.3) 

The value of k associated with the eigenvalue emerging from the continuum solves 
this equation and by Proposition 13 . II it tends to zero. Using these two facts we will 
deduce the asymptotic formula stated in Theorem 12. 21 First of all, in the following 
we will consider the equation ()7.3)1 for real positive n only. This restriction can 
be justified easily, since for negative k the associated function u given by ()4.6p 
increases at infinity and thus it does not belong to L^(n). In order to calculate 
the asymptotics, we extract the leading part of the second term in the equation 
()7.3|1 : for small positive k. we have 

{^P, Ti3(«:, l){I+TsiK)Tn{n, /))"'</)) ^^m) = ^.^in, /)</')^2(n) + T^,{k, I) , (7.4) 
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where Ti3(k, Z) : ^ M is a function defined for (k, /) G x [Aq, +00) which 
satisfies the relation 

T,,M = O + e-'^A (7.5) 

\ cosn kl j 

as (k, /) (0, +CX3). To get this estimate one has to employ the relation Ti4(k, /) = 
C»(||Ti3(K,/)f ) and the fact that 

implied by the definition of T13 - see ()5.9|1 . Our next step is to extract the leading 
term from the first summand at the right hand side of ()7.4|1 . We will do it in the 
same way as in last section, the only difference is that now we have to take into 
account also the second transverse- mode contribution to ()5.9p . 
We introduce the function V\ that is an analogue of ()6.8|) by 

— -Xismx2 / ti smt2 0(i) " ^ , a;i < , 

V^{x) ■= { ^cosh/c/ y_ ^7 7) 

0, xi>0. 

Let W\ be a solution to the problem ()4.5p with v = V\. We also introduce the 
function 'V2 in the following way 



■ sinh \f?)X\ sin 2x2 / sinh a/S^i sin 2t2 0(t) dh, xi < , 



Vgfx) := < TTv^S cosh -\/3/ 

, xi > 

(7.8) 

and suppose that is a solution of ()4.5|) with f = V2. One can check that 
(V^,Ti3(«:,/)0)^2(n) = - (V^,(A + 1) (^V^ + ^(W^-V^)))^^^^^ + T15 («:,/) , 

where V = V1 + V2, W = Wi + VF2, and the function Ti^{k, I) satisfies the estimate 

ri5(«:, l) = (^^^ + + e-2-^'^ (7.9) 

as (k, /) — s> (0, +00). Calculating the inner product {ip, Ti^{k, /)0)j;^2(n) in the same 
way how we deduced ()2.Hj) and bearing in mind the asymptotics 1)2.11) for ip together 
with ()7.4j) . ()7.5j) . and ()7.7j) - ()7.9|) we obtain the equation 



kI 



2K+p—;= — - / tismt2(f)dt 

Jt[ cosh hX 



(7.10) 



+p -j=r- I sinh -\/3ti sin 2^2 4> dH + Tiq{k, /) = 

cosh V 3/ J 



n- 
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where p is again the parity of ip. The function Ti6(k, /) satisfies 

+ + «:^e-2v/^^ + e-^^' (7.11) 

cosh kI cosh KL J 

as (k, /) — i> (0, +00). Since the function obeys ■?/' = T3(/t = 0)0, we can take into 
account the definition of the last operator (see ()4.1H|l ) and the asymptotics ()2.H) 
to conclude that 



p\ — = / ti sin ^2 0(^) d t , 

V vr n J 
n- 

pP = -j= J sinh v^^i sin 2^2 (j){t) d^t . 



which together with ()7.10|) leads us to (/5 = /?„) 



or equivalently, 



2k ; — — + Tiq{k, I) = . 

cosh 2 cosh V 3/ 



■ T^e{n,l) = 0. (7.12) 



cosh 2 cosh ^3/ 



We know that this equation has a positive solution tending to zero as k — ^ 0. In 
view of ()7.11|) . ()7.12|) . and the trivial inequality 

e^ 

1 < — r- < 2 
coshr 

we have for this solution the following estimate, 

<K< C2e-2^' 

with constants Ci, C2 independent on /. Using it we can expand the first term in 
the equation ()7.12|) with respect to nl, which is small, and to estimate T15 - see 
()7.11|) . In this way we arrive at the relation, 

«:-/327rV3e-2^' + 0(e-2^') =0, 

which implies the sought asymptotical expansion ()2.6|) . ()2.7|) . The second formula 
for p stated in Theorem 12.21 can be proven completely by analogy with the proof 
of ()2.5|1 . One just should multiply the equation (A + l)ip"' by e^"^^ sin 2x2 and 
integrate then by parts over 11^ passing then to the limit as i? — +00. 

The argument concerning the asymptotics for the associated eigenfunction is 
completely analogous to that of the previous section. The solution to the equation 
()7.H) is given by 

$ = (/ + T8(k)Ti3(«:, /))->■ 
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Now one has just to perform the expansion using the fact that the operator 
Tg,{K)Ti3{K,l) is small, then using the obtained asymptotics for k, to apply to 
the remainder the estimate ()7.6p . to construct the corresponding eigenfunction of 
the problem ()2.9|1 by the scheme described in the Section 5, and finally, to recover 
the eigenf unctions of Hi{an)- This completes the proof of the second theorem. | 
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